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ABSTRACT 

Perturbative renormalization group theory is developed as a unified tool for 
global asymptotic analysis. With numerous examples, we illustrate its applica- 
tion to ordinary differential equation problems involving multiple scales, bound- 
ary layers with technically difficult asymptotic matching, and WKB analysis. In 
contrast to conventional methods, the renormalization group approach requires 
neither ad hoc assumptions about the structure of perturbation series nor the use 
of asymptotic matching. Our renormalization group approach provides approx- 
imate solutions which are practically superior to those obtained conventionally, 
although the latter can be reproduced, if desired, by appropriate expansion of 
the renormalization group approximant. We show that the renormalization group 
equation may be interpreted as an amplitude equation, and from this point of 
view develop reductive perturbation theory for partial differential equations de- 
scribing spatially-extended systems near bifurcation points, deriving both ampli- 
tude equations and the center manifold. 
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1. Introduction 



Asymptotic and perturbative analysis has played a significant role in applied 
mathematics and theoretical physics. In many cases, regular perturbation meth- 
ods are not applicable, and various singular perturbation techniques must be 
used. Examples of widely-used techniques for ordinary differential equations 
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(ODEs) include ' the methods of multiple scales, boundary layers or asymptotic 
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matching, WKB, stretched coordinates, averaging, the method of reconstitution, 

and center manifold theory. 6 Although these methods are well known, each has 
its own drawbacks, preventing mechanical (or algorithmic) application. Indeed, 
it is probably fair to say that the practice of asymptotic analysis is something of 
an art. 

Multiple- scales analysis has proven to be a particularly useful tool for con- 
structing uniform or global approximate solutions for both small and large values 
of independent variables. In this method a set of scaled variables, which are re- 
garded as independent variables (although they are ultimately related to one 
another), is introduced to remove all secular terms. The choice of the set is, in 
some cases, nontrivial, and may only be justified post hoc. Nevertheless, this 
method is usually considered the most general, subsuming the others mentioned 
below. 

Differential equations whose highest order derivatives are multiplied by a 
small parameter e often yield solutions with narrow regions of rapid variation, 
known as boundary layers. Boundary layer techniques can be applied if the 
thickness of these regions tends to zero as e — > 0; otherwise, WKB must be used. 
The limitation of WKB is that it applies to linear equations only. Although 
boundary layer methods apply to nonlinear as well as to linear problems, the 
determination of the expansion parameter can be subtle. Furthermore, matching 
of outer and inner expansions via intermediate expansions is required, sometimes 
involving delicate arguments that are difficult to perform mechanically. 

Another class of related problems concerns partial differential equations 
(PDEs) describing nonequilibrium, spatially-extended systems near bifurcation 
points. Such systems often exhibit spatial-temporal patterns modulated by an 
envelope function (or amplitude) which varies slowly compared with the pattern 
itself. Extracting the long wavelength, slow timescale behavior of such systems 
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is the task of reductive perturbation methods, which are themselves related to 
multiple-scales analysis. 

The purpose of this paper is to present a unified, and physically-motivated ap- 
proach to these classes of problems, based upon the renormalization group (RG) . 
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The essence of the renormalization group method is to extract structurally sta- 
ble features of a system which are insensitive to details. For example, field 
theories, critical phenomena, polymers and other statistical mechanical systems 
exhibit universal scaling functions and critical exponents in the limit A/£ — > 0, 
where A is some ultra-violet cut-off and £ is the (temperature-dependent) cor- 
relation length. The renormalization group is the principal tool with which to 
elucidate this universal behavior and is properly regarded as a means of asymp- 
totic analysis. 
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The usefulness of this point of view has been amply demonstrated by the 
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relationship between the renormalization group and intermediate asymptotics. 
In particular, the large-time asymptotic behavior of certain initial-value prob- 
lems is given by a similarity solution of the governing PDE, where the similarity 
variable contains anomalous exponents which may not be determined a priori 
by elementary dimensional considerations. Nevertheless, renormalized perturba- 
tion theory combined with the renormalization group, gives an expansion for the 
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anomalous exponents and the solution. 

The similarities between renormalization group and singular perturbation 
methods extend also to technical details: both perturbative renormalization 
group and conventional singular perturbation methods remove secular or diver- 
gent terms from perturbation series. These formal similarities invite a natural 
question: what is the relation, if any, between conventional asymptotic methods 
and the renormalization group? 

In this paper, which is an extended version of our preliminary report, 16 we 
demonstrate that singular perturbation methods may be naturally understood 
as renormalized perturbation theory, and that amplitude equations obtainable 
by reductive perturbation methods may be derived as renormalization group 
equations. 

Our studies indicate that the renormalization group method may have several 
advantages compared with conventional methods. Although we recognize that 
our analysis is at the formal, heuristic level, we suggest that a more careful 
mathematical analysis would be worthwhile, given the potential usefulness of our 
central claim. 

One advantage of the renormalization group method is that the starting point 
is a straightforward naive perturbation expansion, for which very little a priori 
knowledge is required. That is, one does not need to guess or otherwise introduce 
unexpected fractional power laws or logarithmic functions of e in an ad hoc man- 
ner. It seems that these e-dependent space/time scales arise naturally during the 
analysis. 
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We will show that the renormalization group approach sometimes seems to 
be more efficient and accurate in practice than standard methods in extracting 
global information from the perturbation expansion. Standard methods often 
attempt to represent an asymptotic solution in terms of asymptotic sequences of 
a few simple functions of the expansion parameter, such as exp, log, powers and so 
on. The renormalization group can generate its own problem-adapted asymptotic 
sequence without matching: in the examples given in section 4, these turn out to 
be complicated functions conveniently defined by an integral representation. For 
small e, this asymptotic sequence can be expanded to reproduce the solutions 
conventionally obtained by asymptotic matching, although in the examples that 
we have studied so far, the conventional approximant is practically inferior to 
the one obtained by the RG. In switchback problems, the RG perturbation series 
may need to be carried out to higher than lowest order, then expanded in e, in 
order to reproduce the (inferior) conventional result. 

A related advantage of the renormalization group seems to be the lack of 
necessity to perform asymptotic matching. To illustrate this assertion, in section 
3 we solve several ODEs with boundary layers, and in section 4 we address the 
difficult technical problem of switchback terms. 

The renormalization group methods for partial differential equations such 
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as the Barenblatt equation, ' ' and front propagation problems in reaction- 
diffusion equations, 11 are, in retrospect, examples of the general approach dis- 
cussed in this paper. We emphasize that our renormalization group method has 
no connection with the so-called method of renormalization or uniformization 
in the conventional perturbation literature; the latter is a mere variant of the 
method of stretched coordinates, and of narrow limited use. 

Lastly, we wish to point out that recently, a method utilizing an invariance 
condition in the solution of multiple- scale singular perturbation problems was 
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proposed independently by Woodruff, based on ideas related to the renormal- 
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ization group. In addition, Kunihiro has shown that the renormalization group 
method that we proposed in ref. 16 may be interpreted geometrically in terms 
of the classical theory of envelopes. 

The outline of this paper is as follows. In Section 2, we discuss the general 
relation between multiple- scale analysis and renormalization group. In Section 
3, we show how boundary layer and WKB problems can be solved using the 
renormalization group. In Section 4, we demonstrate with several examples that 
the renormalization group approach has technical advantages to conventional 
asymptotic methods. In Section 5, the renormalization group is applied as a 
reductive perturbation tool to the derivation of global slow motion equations for 
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partial differential equations. Center manifold theory is also briefly considered 
from the same point of view. We conclude in Section 6. 



2. Multiple Scale Theory and RG 



In this section, we show that multiple-scale analysis is equivalent to the RG, 
and that the solvability condition used in multiple scales to remove the secular 
divergences is equivalent to the physical assumption of renormalizability in RG 
theory. 

2.1. Rayleigh Equation 
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The example we consider below is the Rayleigh equation, closely related to 
the Van der Pol oscillator: 



It is known that the method of uniformization or renormalization fails here, 
and this example is a textbook illustration of multiple scales analysis. We show 
here that from only the simple-minded straightforward expansion, not only is the 
RG capable of identifying automatically all different multiple scales required by 
multiple scales analysis, but also produces a uniformly valid asymptotic solution 
without encountering the ambiguity which often plagues higher order calculations 
in multiple scales analysis. 

A naive expansion y = yo + ey\ + e 2 y2 + ■ ■ • gives 



where Rq, Go are constants determined by the initial conditions at arbitrary t = 
to. This naive perturbation theory breaks down when e(t — to) > 1 because of 
the secular terms. The arbitrary time to may be interpreted as the (logarithm of 
the) ultraviolet cutoff in the usual field theory. 11 To regularize the perturbation 
series, we introduce an arbitrary time r, split t — to as t — r + r — to, and absorb 
the terms containing t — to into the renormalized counterparts R and G of Ro 
and Go, respectively. This is allowed because Ro and Go are no longer constants 
of motion in the presence of the nonlinear perturbation. 




(2.1) 



l 
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(2.2) 
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We introduce a multiplicative renormalization constant Z\ = 1 + a n e n 
and an additive one Zi = bn£ n such that Ro(to) = Zi(to,r)R(r) and 
©o(^o) = ©( T ) + Z2(to,r). The coefficients a n and 6 n (n > 1) are chosen or- 
der by order in e to eliminate the terms containing r — to as in the standard 
RG. The choice a\ = -(1/2) (1 - R 2 /4)(r -t ),bi = removes the secular 
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terms to order e, and we obtain the following renormalized perturbation result 

y(t) ={R + e| (l - ?pj (t-r)} sin(t + 0) 

1 R 3 ^ 

- e—R 3 cos(t + 6) + e— cos 3(t + ©) + 0(e 2 ), 

where i?, are now functions of r. Since r does not appear in the original 
problem, the solution should not depend on r. Therefore, {dy/dr)t = for 
any t. This is the RG equation, which in this case consists of two independent 
equations 

dR 1 „ A 1 „<A . 2 x dO _ 2 



_ =^1-^+0^), _ =OM . ,2.4) 
Solving (2.4) and equating r and t eliminates the secular term, we get 

R(t) = R(0)/^e-^ + ^(0) 2 (1 - e- rf ) + 0(e 2 t), 6(t) = 0(0) + 0(e 2 t), (2.5) 

where i?(0),©(0) are constants to be determined by the initial condition. As- 
suming the initial condition y(0) = 0,y'(0) = 2a, we find R(0) — 2a, ©(0) = 0, 
and the final uniformly valid result reads 

y(t) = R(t) sin(t) + |r^) 3 {cos(3t) - cos(t)} + 0(e 2 ), (2.6) 

which approaches a limit circle of radius 2 as t — > oo. 

The second order RG calculation shows the assumption of perturbative renor- 
malizability is consistent and no ambiguity arises at all. The corresponding am- 
plitude and phase equation to order 0(e 3 ) are 



/ R*\ (2 - 7) 

l-Sr)+0(e 3 ), 



dt 2 

d© _ e 2 / R 
~dt ~ ~ ~8 V ~ 32 , 

from which the multiple time scales T\ = e£, T2 = e 2 t, • • • used in multiple scales 
analysis appear naturally (although the RG does not require such identifications) . 
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When R = 2, (2.7) reduces to 



dR 
~dt 



+ O(e 3 ), 



~dt 
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e 2 + 0(e 3 ). 



(2.8) 



In this simple example, it was straightforward to determine the multiple time 
scales. However, it is well known that in many cases, within multiple scales anal- 
ysis hidden intermediate scales must be included in the perturbation expansion 
so as to obtain the correct result. In the next example, will show that the RG 
method is a more straightforward but secure way to determine multiple slow time 
scales than the multiple scales method. 

2.2. Mathieu Equation 

The second illustrative example we examine using RG is the Mathieu equa- 



where a and e are parameters. 

The Floquet theory of linear periodic differential equations 1 predicts that 
in the (a, e) plane there are some regions where the solutions to (2.9) remain 
bounded for all t and stable, and others where the solutions are unstable. Per- 
turbative investigation shows that for sufficiently small e, all solutions y(t) are 
stable for a > 0,a / n 2 /4, n = 0, 1,2, Without loss of generality, we in- 
vestigate the stability of solutions near a = 1/4 and e = to find the stability 
boundary in the (a, e) plane. We treat the boundary curve a as a function of e 
and expand a in powers of e: a(e) = l/4 + aie + a2e 2 + • • •• It is our goal to deter- 
mine values of ai, a2, • • • perturbatively. Multiple-scale analysis can be applied to 
this problem, and the coefficients a\ = l,a2 = —1/2 are determined. However, it 
turns out that the introduction of multiple time scales r\ = et, T2 = e 2 t, • • • is not 
sufficient to determine the second order coefficient ai even after the first order 
coefficient a\ is set to 1. Through careful analysis, it is found that a new hidden 
time scale a = e 3 / 2 t must be introduced into the problem, and the perturbative 
expansion must be done in powers of e 1 / 2 , rather than the original expansion in 
powers of e. It is necessary to go to the fourth order in powers of e 1//2 to determine 
<22- Thus, the procedure required to determine all necessary time scales is not 
mechanical: if any hidden scales are omitted or cannot be determined, correct 
results will not be guaranteed. This represents a typical shortcoming of multiple 
scales analysis. 



tion 
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+ (a + 2ecost)y = 0, 



(2.9) 



dt 2 
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Now we demonstrate how the unexpected time scales such as a = e>lH ap- 
pear automatically from the RG equation, starting only with a straightforward 
perturbative expansion. Substituting a = 1/4 + a\e + ci2e 2 + ■ ■ ■ in (2.9) and 
expanding in powers of e (not e 1 / 2 ) as y = yo + ey\ + e 2 y2 + ■ ■ •, we get 

~dW + 4 m = °' (2 - 10) 
^ + l yi = -(ai + 2 cos*)y , (2.11) 

^^ + ^y2 = -a2yo-(ai + 2cost)yi, (2.12) 

and so on. First, let us determine the first order coefficient a\. The straightfor- 
ward perturbation result, to 0(e), is given by 

y(t) =R cos(t/2 + 6 ) +eR \ -I cos(t/2 + O ) + I cos(3t/2 + 6 ) 

1 2 2 (2.13) 

- ai(t - t ) sin(t/2 + 6 ) - (t - to) sin(£/2 - O )} + 0(e 2 ), 

where Ro, Oo are constants dependent on initial conditions given at some arbi- 
trary time to- Similarly, the secular divergences can be removed by regarding 
to as a regularization parameter and renormalizing the bare amplitude Aq and 
bare phase Go: -Ro(^o) = Zi(to, (j)R((~i), ©o(^o) = -^2(^0? A 4 ) + ©(/•*)) where \i is 
some arbitrary time scale, as was done in previous problems. The renormalized 
perturbation result is 

y(t) ={r(/u) + eR( - 1/2 + (t - fj,) sin 20^)) ) cos(t/2 + 0) - eR 

R o ( 2 - 14 ) 

x (ai + cos20)(t- / u)sin(t/2 + 0) + e-cos(3t/2 + 0) +0(e 2 ). 

The RG equation dy/dji = for any t gives 

= e#sin20 + O(e 2 ), ^ = e(a x + cos 20) + 0(e 2 ). (2.15) 

For convenience, we introduce the complex amplitude A = Re lQ as A = B + iC, 
with its real and imaginary parts B = i?cos6, C = i?sin6. The equations for 
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B(n) and C(n) are 



e(l - ai)C(//), C"(//) = e(l + ai)B(//). 



(2.16) 



Thus, we have 



B"^) = e 2 (l-a\)B^). 



(2.17) 



Solving this and setting ji = t, we get 



B(t) = K x e ± ^~< 



(2.18) 



where K\ is a constant, and the first slow time scale r\ = et has appeared 
automatically. Obviously, for |ai| < 1, instability sets in, where the solution 
grows exponentially with time t, while for |ai| > 1, the solutions are bounded 
and stable. Therefore, near e = 0, the stability boundary isa = l/4±e + 0(e 2 ). 

We now set a\ = 1 and go to the second order to determine 02- For order e 2 , 
a special solution to (2.12) is obtained 



Extending the renormalization procedure to the second order, we find all the 
secular divergences to this order can be removed completely, a sign of the consis- 
tency of perturbative renormalizability. Keeping only the two lowest harmonics 
with prime frequency and omitting other higher frequency terms which are not 
important for determining the stability boundary, we obtain the renormalized 
perturbation result, to order e 2 , 



where H.F.T. represents all higher frequency terms. The RG equation to order 




(2.19) 




cos 20) 



(2.20) 
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e 2 now reads 

^ = e J Rsin2e + 0(e 3 ), ^ = e(l + cos 20) + 1/2) + 0(e 3 ). (2.21) 

Accordingly, the equations for £?(//) and C(/u) become 



B'M = - 
Thus, we get 



e 2 (a 2 + l/2) 
1 - e/2 



2e + 



e 2 (a 2 + 1/2) 
1 - e/2 



(2.22) 



2e 3 (a 2 + l/2) e 4 (a 2 + 1/2) 



1-6/2 (l-e/2) 2 
Keeping only the lowest order term of (2.23) gives 

B"{n) w -e 3 (2a 2 + + 0(e 4 ) 

which has the solution (setting fj, = t) 

B(t) = K 2 e ± ^ Tl " 3/2 \ 



(2.23) 



(2.24) 



(2.25) 



where K 2 is a constant, and the second and the third slow time scales a = 
e 3//2 t, r 2 = e 2 t appear naturally. We apparently have stable solutions for a 2 > 
— 1/2 and unstable solutions for a 2 < —1/2. Therefore, to order e 2 , the instability 
boundary is given by 



o(e) 



^ + e-^e 2 + 0(e 3 ), e 



0. 



(2.26) 



2.3. Oscillator with Time Dependent Spring Constant 

28 

The third illustrative example is an oscillator governed by the equation 

^ + y-ety = 0. (2.27) 

The initial conditions are y(0) = 1 and y'(0) = 0. The regular perturbation 
theory breaks down for t — > oo, and multiple scales analysis can be applied to 
eliminate the secular behavior. However, it turns out that multiple time scales 
must be chosen as to = t, t\ = e 1//2 t, r 2 = et, ■ ■ •. Since the frequency of the 
oscillator is found to be time-dependent, the method of stretched coordinates or 
the so-called method of uniformization or renormalization (in the conventional 
applied mathematics sense) does not work here. 
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We will see how a uniformly valid solution can be constructed simply from 
the naive perturbation series with the aid of the RG. To solve (2.27), we assume 
a straightforward expansion in powers of e (not e 1 / 2 ), y = yo + ey\ + e 2 y2 + • • •• 
The bare perturbation result, to order e, is given by 

y(t) = R cos(t + Go) + eRofyt 2 - t 2 ) + ^(t - t )} sm(t + G ) + 0(e 2 ). (2.28) 

As in the preceding examples, renormalizing the bare amplitude Rq and phase 
Go removes the secular divergences. The renormalized perturbation result is 

y (t) = | J R+^e J R(t-^+ai)}cos(t+G)+^e J R(t 2 -^ 2 +6i)sin(t+G)+0(e 2 ), (2.29) 

where R, G are functions of arbitrary time scale fx, and a\, b\ are arbitrary con- 
stants. The RG equation reads 

Solving (2.30) and setting \i = t in (2.29) give 

R(t) = R(0)e^ et + 0(e 2 t), Q(t) = - ^et 2 + 0(0) + 0(e 2 t). (2.31) 

Thus, we obtain the uniformly valid result 

y(t) = R(t)cos(t + e(t)) + -eR(t) (ai cos(t + G) + h sm(t + G)) +0(e 2 ). (2.32) 

Imposing the boundary conditions y(0) = l,y'(0) = gives -R(O) = 1,0(0) = 
0, a\ = 0, bi = — 1. Therefore, the final result is 

y (t) = e ^ et cos(t - jtt 2 ) - ^ee3 rf sin(t - ^et 2 ) + 0(e 2 ), (2.33) 

where the frequency defined as uj = dQ/dt becomes time-dependent: uj = 1 — 
^et + 0(e 2 ). Rewriting et 2 as (e 1//2 t) 2 , two slow time scales T\ = e 1//2 t, T2 = et are 
easily identified from the RG result (but these identifications are unnecessary in 
our approach). 

The RG scheme given above is also applicable to quantum systems with 
discrete or continuous energy spectrums, especially those which involve resonance 
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phenomena, e.g., the Rabi flopping, the Stark shift, the Bloch-Siegert shift. 
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The multiple time scale perturbation analysis has successfully given a unified 

30 

framework for all quantum resonance. In a similar way, the RG method simply 
recovers all resonance equations which turn out to be simply RG equations. The 
application of RG to the time-dependent Schrodinger equation also reproduces 
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the Fermi's Golden Rule. Here we will not give detailed calculations of these 
problems. In the next section, we will show that WKB problems can be easily 
solved using the RG method. Therefore, many quantum problems which are 
usually solved using WKB and/or multiple scales analysis can also be studied 
using the RG approach. 

To summarize, it seems that the RG method is more efficient and mechanical 
than the multiple scales method in determining the multiple slow time scales. In 
the RG approach, the starting point is simply a straightforward naive perturba- 
tion series, and all necessary multiple scales arise naturally from RG equations. 
The above examples reveal two important points, demonstrated more generally 
below: (1) the results of multiple scales analysis can be obtained from renormal- 
ized perturbation theory, and (2) the RG equation describes the long time scale 
motion of the amplitude and the phase. 

3. Boundary-Layer Theory, WKB and RG 

Another important class of singular problems is that for which the highest 
order derivative of the equation is multiplied by a small parameter e, e.g., WKB 
and boundary layer problems. 

Boundary-layer theory and asymptotic matching are a collection of singular 
perturbation methods for constructing a uniformly and globally valid solution 
by calculating the separated outer and inner solutions and then matching them 
across intermediate scale solutions. Quite often, the intermediate matching is 
very lengthy and only some particular matching method will work. WKB theory 
is well known to be a powerful tool for obtaining a global approximation to 
solutions of a linear differential equation whose highest derivative is multiplied 
by a small parameter e. Many linear problems often solved by WKB theory can 
be solved by boundary layer theory; indeed, in these cases, boundary layer theory 
(thickness of the boundary layer goes to zero as e — > 0) is a special case of WKB 
(thickness of the boundary layer remains finite even as e — ■> 0) . The limitation of 
the conventional WKB method is that it applies only to linear problems, while 
boundary layer theory works for linear as well as nonlinear problems. 

In this section we will demonstrate explicitly that many boundary layer 
problems, linear or nonlinear, can be solved by the RG. The uniformly valid 
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asymptotics of boundary layer problems can actually be constructed from the 
inner expansion alone, with the aid of the RG, without the need for intermediate 
matching. 



3.1. Simple Linear Example 



Consider the following simple example, which describes the motion of an 
over damped linear oscillator: 

«|? + * + »-0, (3.1) 

where e is a small parameter. A standard dominant-balance argument shows that 
there exists a boundary layer of thickness 5 = 0(e) at t = 0. Thus, we set t = er, 
and rewrite equation (3.1) as 

£ + * + *-a (3,, 

Naive expansion gives 

y(r) = Aq + B e~ T + e[ - Aq(t - tq) + Bq(t - r )e- T ] + 0(e), (3.3) 

where the coefficients Aq, Bq are constants of integration and 0(e) refers to 
all the regular terms of order e and higher, which are finite even in the limit 
t — tq — > oo. This naive perturbation theory breaks down due to the diver- 
gence of secular terms for large r — tq. However, this divergence can be re- 
moved by regarding tq as a regularization parameter and renormalizing Aq,Bq 
as Aq(tq) = Z\A(\i), and Bq(tq) = ZiB([i). Here fi is an arbitrary time, and 
A, B are the renormalized counterparts of Aq, Bq. The renormalization constants 
= Etf ' a n(rQ, ^)e n , Z 2 = Ya? b n (TQ, n)e n (a = l,b = 1) are chosen order by 
order in e to eliminate the secular divergences. Split r — tq as (r — /i) + (fx — tq), 
and then absorb the divergent part fx — tq in the limit tq — > — oo by redefining 
Aq and So- Choosing a\ = il — TQ,b\ = —(fx — tq), we get the renormalized 
perturbation result 

y(r) = A(jjl) - eA(fi)(r - /*) + [B(fu) + eB(fu)( T - v)]e~ T + 0(e). (3.4) 

However, it is impossible that the actual solution y(r) can depend on the ar- 
bitrary time [i which is not present in the original problem. Thus we have the 



13 



renormalization group equation dy/dfx = for any r, which gives 



dA 

— + eA + 
djJL 



d[i 



e- r + 0(e 2 ) =0, 



or 



dA n dB 9 

dfi dfi 



(3.5) 



(3.6) 



Extending the RG calculation to the second order gives, without any ambiguity, 
-(eA + e 2 A) +0(e 3 ), ^ = eB + e 2 B + 0(e 3 ). (3.7) 



cL4 



Solving them, setting fx = r and setting back r = i/e in (3.4), we finally obtain 
the uniformly valid solution 



y (f) = Cie-( 1+ ^ + C 2 e- t / e+ ( 1+£ ) t + 0(e 2 ), 



(3.8) 



where C\,Ci are constants to be determined by the initial conditions. Clearly, 
the RG result to order e 2 recovers exactly that obtained by the standard singular 
methods. Notice that the equations in (3.7) are nothing but the equations of 
motion for slow time scale: the amplitude equations. Thus, amplitude equations 
are renormalization group equations. We announced this result previously, and 
derived the Burgers equation as a renormalization group equation. A much 
more complicated example illustrating this point will be given in Section 5. 



3.2. Example with log e 

The second example we consider is 32 

ey" + xy' - xy = 0, y(0) = 0, y(l) = e. (3.9) 

A standard dominant-balance argument tells us that there exists a boundary 
layer of thickness of order e 1 / 2 (but not e) at x = 0. The complication in the con- 
ventional asymptotic matching stems from the fact that the inner expansion must 
contain not only powers of e 1//2 but also those terms containing combinations of e 
and log e to make the intermediate matching successful. Here we explicitly show 
that the renormalized naive inner expansion in powers of e 1//2 gives a uniformly 
valid asymptotic solution. This reveals that those unexpected terms containing 
log e in the conventional approach are just an artifact of perturbative expansions 
of x~ e . 
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Assuming x = e 1//2 X, and y(x) = Y(X), we transform (3.9) into 
d 2 Y dY 

IX^ +X IX~ VlXY = ^ nO) = 0,Y(lM) = e. (3.10) 
Naive expansion in e 1 / 2 , Y(X) = Y (X) + e^Y^X) + eY 2 (X) + ■■■ gives 

r " + Xy ' = 0, Y^ + XY' n = XY n _ x , (n>l), (3.11) 

Thus the naive perturbation result to order e is 

x x 
Y(X) ~A + B J dse-* 2 ' 2 + € 1 / 2 {Aq{X -X Q ) + B (X - Xq) J dse^ 2 / 2 



X 

+ R.T.} + e{^A (X - X ) 2 + ^B (X - X ) 2 J dse-* 2 ' 2 

o 

X 

-(^» +B »Mf)/^ /2+R 4 



(3.12) 

where Aq,Bq are integration constants, and R.T. represents regular terms finite 
even in the limit X — Xq — > oo and \og(X / Xq) — > oo. The divergence can 
be controlled by renormalizing Aq = ZiA(p), Bq = Z2B(p>), where Z\(p) = 
a n e n / 2 , ao = 1 and Zi{\i) = X]o° b n £ n ^ 2 : , bo = 1, are renormalization constants 
and fi is some arbitrary position. The choice a\ = Xq — fj,, ci2 = (1/2) (Xq — (i) 2 
and b\ = Xq — ^,62 = {1/2)(Xq — /u) 2 — (-^A + B) log(Xo/ fj,) successfully removes 
the divergences up to order e, and the renormalized perturbation result is 

Y(X) ~{^) + e x l 2 A(X - ii) + t X -A(X - ^) 2 } + [b^i) + e l ' 2 B{X - fi) 

x 

+ e^B(X - /j) 2 - e(^=A + B)log(X//j)} J dse^' 2 . 



(3.13) 

e 1/2 A + 0(e 3/2 ), (3.14) 



The RG equation dY /dji = gives 

dA _ 



d 4- = e 1/2 S - e(^=A + + 0(e 3 / 2 ). (3.15) 
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Solving these two equations, we obtain 

A(fx) =Cie £l/ ^ + 0(e 3 /V), 



(3.16) 



B(ji) = ~ ^C^^ + C 2 yr e e^» + 0(e 3 /V), (3.17) 

where Ci, C2 are constants to be determined by the given boundary conditions. 
Setting /j, = X, we obtain 

x 

Y(X) ~ C ie £l/2x + { - ^ + C 2 X- e }e £l/2x J dse~ s */\ (3.18) 



Imposing boundary conditions Y(0) = 0,Y(l/y/e) = e gives C\ = and C2 = 
V^TTre -6 / 2 = v /277re^ eloge as e -> 0+. Setting back X = x/e 1/2 , we obtain the 
final uniformly valid asymptotic result to order e, 

00 

y{x) ~ e^- £ [l - y dse~ s2/2 }. (3.19) 

Thus the terms such as e log e, which are present in the inner expansion given in 
{e.g.,) ref. 1, are relics of the expansion of x~ e . It is worthwhile to note that 

the RG result is slightly different from the asymptotic matching result given by 

32 

Bender and Orszag in their book. To leading order, the former is 

00 

y§ G {x)~e x {\-y/2fi J rfse^ 2 / 2 }, (3.20) 

while the latter is 

00 

yg°(x) ~ e x - J dse-* 2 ' 2 . (3.21) 

Comparing with the numerical result of the original equation (3.9), we find that 
in the boundary-layer region, the RG result (3.20) is a better approximant than 
the standard result (3.21). 
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3.3. Nonlinear Boundary Layer Problem 



Boundary-layer analysis applies to nonlinear as well as to linear differential 
equations. In this section and in the following section, we will demonstrate that 
the RG method can be used to solve nonlinear boundary layer problems. 

33 

Let us consider the following illustrative nonlinear problem: 

ey" + 2y' + e» = 0, y(0) = y(l) = 0. (3.22) 

There is only one boundary layer of thickness e at x = 0. Setting X = 
x/e,Y(X) = y(x) in (3.22) gives 

d 2 Y dY Y . . 

«i + 2 lx = - < X23) 

Assuming an inner expansion Y = Yq + eY\ + ■ ■ ■ gives the following asymptotic 
result as X — > oo, 

F(X) ~ A + S e- 2X - e!^e Ao {X - X ) + R.T. J + 0(e 2 ), (3.24) 

where Aq, Bo are integration constants and R.T. represents all regular terms in the 
expansion finite even in the limit X — X$ — > oo. The renormalized perturbation 
result obtained as in the previous examples is 

Y(X) ~ A(p) + B(p)e- 2X - e^e A M(X - fi) + 0(e 2 ). (3.25) 

The RG equation gives, to order e, 

dA 1 A dB . 

Solving (3.26) , we get 

where Ci, C2 are constants of integration to be determined by the given boundary 
conditions. Equating fx and X in (3.25) and restoring x = eX, we obtain the 
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uniformly valid asymptotic result 



y(x) ~ log (jf^j + C 2 e- 2x ^ + 0(e). (3.28) 

Imposing boundary conditions y(0) = 0, 2/(1) = gives C\ = 1,6*2 = - log 2 in 
the limit e — > 0+. Therefore, the final result is 

y(x) ~ log (^|y) " (log2)e- 2a; / e + 0(e). (3.29) 
This RG result recovers the leading order result from boundary layer analysis. 



3.4. Nonlinear Problem of Carrier 

34 

In this section, we consider a first-order nonlinear model problem of Carrier, 
(x + e/)/' + f=l, f(l) = 2, 0<x<l. (3.30) 

The exact solution can be obtained by integrating (3.30) once, 

fM = -f + g + ?fe±il + 4) . (3.31) 

It becomes, however, a nontrivial singular perturbation problem, if we pretend 
that we cannot obtain the exact solution. The method of strained coordinates 
or the method of asymptotic matching can be applied with a rather lengthy 
matching. We show here how to solve the problem using RG without matching, 
and give the exact result, starting only from the inner expansion. 

First, we apply the usual dominant balance argument to make the structure 
of the equation clear. We introduce X = 77(e) a; and F = 8(e) f; the latter is 
needed because the equation is nonlinear. The original equation reads (X + 
(erj/S)F)dF/dX + F = 5. erj/S 1 corresponds to the outer limit, and ei]/5 ~ 1 
is the only nontrivial alternative possibility. Hence, 5 ~ e a and 77 ~ e a -! with 
a G (0, 1] are the useful scalings. The expansion parameter becomes 5 ~ e a . It 
turns out that any choice of a is admissible in this case, and so we adopt the 
simplest choice a = 1. 



18 



Accordingly, we rescale / as F = ef to convert the original equation (3.30) 

to 

(x + F)F' + F = e. (3.32) 
Expanding F as F = Fq + eF\ + • • •, we have 

(a; + F Q )F' Q + F = 0, (3.33) 

whose general positive solution is 

F Q (x) = (x 2 + Aq) 1 ' 2 - x, (3.34) 

with Aq a constant of integration determined by the initial condition given at 
some arbitrary xq. The first order equation is given by 

(x + F )F[ + F\Fq + F - 1 = 1. (3.35) 

This linear equation has a general solution 

F ^ = 1 Z'aW (3 - 36) 
[x l + Aq) l I z 

Thus, the straightforward perturbation result, to 0(e), is given by 



F(x) = (x 2 + A Q fl 2 - x + e [ X ; o) 2 + 0(e 2 ). 

[X l + Aq) 1 / 1 



(3.37) 



We see that this naive perturbation (3.37) breaks down formally for x ^> xq. 
Actually, the domain of our problem is finite, and because x is not scaled, it is 
not possible that x ^> xq can occur within the domain. A better argument is as 
follows. Since the boundary condition is F(l) = 2e, near x = 1 the 0(e) term 
dominates; this is a singular perturbation, and indeed the perturbation term 
diverges relative to the zeroth order term. 

The secular divergence can be removed by renormalizing Aq by Aq(xq) = 
ZA(fi), and the renormalized perturbation result obtained is 

F(x) = (x 2 + A(^' 2 - x + * { Jl~£ ))1/2 + 0(e 2 ). (3.38) 

The RG equation gives, to 0(e), 

dA/d/j = 2e (3.39) 
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with solution 

A(p) = A(0) + 2e/j. (3.40) 
Setting /j, = x and / = F/e, we obtain the uniformly valid asymptotic result 

,/ n x fx 2 2x A(0)\ 1/2 . . 

/ ^ = ~ + U + T + -? i J • (3 - 41) 

Imposing the boundary condition /(l) = 2 gives A(0) = 2e + 4e 2 . Therefore, the 
uniformly valid result to order e, is given by 

*t \ x ( x2 2(x + l) A 1/2 ,„ „ nN 

/(x, e) = --+ ( pr + 1 - J +4j . (3.42) 

This happens to be the exact solution to the problem. A further calculation 
demonstrates that all the higher order corrections vanish. The conventional 
methods can also recover the exact result, but clearly the RG is simpler. 

3.5. Problem with Multiple Boundary Layers 



In many situations, there exist multiple boundary layers at one side, for 
which multiple calculations of inner and outer solutions and their asymptotic 
matchings have to be made in different separated regions to obtain a uniformly 
valid solution. Again it turns out that the RG method manages to produce the 
solution without any matching needed. Let us consider the following initial-value 
problem 

e vy» + ( e i/2 

+ e + e 3 / 2 )y" + (1 + e 1/2 + e)y' + y = 0, (3.43) 

with initial conditions y(0) = 3,2/(0) = -1 - e~ 1/2 - e -1 ,y"(0) = 1 + e" 1 + 
e -2 . The exact solution is y(t) = e~ x + e~ x ^ l/2 + e~ x l e . Pretending we do not 
know how to solve it exactly, we resort to conventional singular perturbation 
methods. It turns out that the conventional perturbation calculation is very 
tedious and rather challenging. By dominant balance, this problem is found to 
have two distinguished boundary layers at t = 0, of thickness of order e 1 / 2 and 
e, respectively. Therefore, one outer solution and two inner solutions must be 
calculated and two asymptotic matchings are necessary, if boundary layer theory 
is used. Starting only with the thinnest or innermost boundary layer by rescaling 
t by t = eT, and expanding y = Y{T) (e.g.,) in e 1 / 2 , the RG method successfully 
recovers the exact solution without any matching. 
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3.6. Linear Boundary-layer and WKB problems I: no turning points 



To conclude this section, we show how linear boundary layer and WKB prob- 
lems in general forms can be treated using RG in a unified fashion. This rela- 
tionship between boundary layer theory and WKB is explained in Bender and 
Orszag's book. The boundary-layer type problem we wish to study using RG has 
the following general form: 

e2 ^ + a ^ X ^~ b ^ y = °' °- a; - 1 ' e ^°+' ( 3 - 44 ) 

where we assume that a(x) is differentiable and b(x) is an arbitrary, not neces- 
sarily continuous function. This equation covers all linear examples we presented 
earlier in this section. A simple dominant-balance argument determines that in 
general, the boundary layer lies at x = when a(x) > for < x < 1, and that 
the boundary layer lies at x = 1 when a(x) < for < x < 1. Without loss of 
generality we will consider only the former case. 

Although in a number of cases we could perform perturbative RG analysis on 
the original general equation (3.44), often it is wiser to start with the canonical 
form of equation (3.44) under the transformation 



y(x) = exp 

converting (3.44) to 



X 

J a{x')dx' 



u(x), (3.45) 



<?^ = Q(x)u(x), (3.46) 

with 

Q(x) = -^a 2 (x) + ±a'(x) + b(x). (3.47) 

This form is just the Schrodinger form, which can be solved by WKB methods. 
Consequently, we can treat both linear boundary layer and WKB problems in a 
unified way. 

In the remainder of this section and in the following section, we will show 
how to solve Schrodinger equations using RG. Our strategy is to first introduce 
a natural change of the independent variable which allows one to obtain effi- 
ciently the non-perturbative part of the solution. The transformation is identical 
to the independent variable portion of the standard Liouville-Green transforma- 
tion or its natural generalization used by Langer, but the crucial difference 
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is that we do not introduce the new dependent variable. This is the analogue 
of the geometrical-optics approximation in WKB theory, and is the starting 
point of a renormalized perturbation series, which reproduces the physical-optics 
and higher-order WKB approximations. Although it may be possible to derive 
even the geometric-optics approximation using RG, we have not succeeded in 
so doing, transformation is sensitive function Q(x). When Q(x) vanishes, its 
zeros lead to turning points in the standard WKB approach. The simplest WKB 
approximations break down there, and connection formulae are required in the 
conventional procedure in order to match approximations on either side, other 

37 

hand, our procedure leads to a uniformly valid approximation. Langer found 
that a suitable generalization of the Liouville- Green transformation can produce 
a uniformly valid approximation across the turning point. Again for the cases 
with turning points, we transform the independent variable only with a straight- 
forward generalization of the no-turning point case. We emphasize that we are 
able to avoid the need to perform matching, and that the transformation of the 
dependent variable is produced naturally by RG. The use of RG is not responsible 
for the choice of the transformation of the independent variable, but our choice 
not to introduce the transformation of dependent variables in contrast to the ap- 
proaches by Liouville and Green and Langer is motivated by RG. This allows us 
to choose a better transformation of the dependent variable, which agrees with 
the conventional result in the small e limit. The corrections and prefactors which 
accompany the zeroth order Langer-type solution are calculated by RG, and do 
differ from and improve upon those obtained by the standard analysis. Further- 
more, we can expand our asymptotic sequence in e to reproduce the standard 
textbook results. 

The remainder of this section concerns Schrodinger problems with no turning 
points. The following section discusses the case with one turning point, and gives 
an outline of how the methods can be generalized for higher numbers of turning 
points and for multiple-boundary-layer linear problems as well. 

We will first rederive the well-known physical-optics approximation using 
the RG theory, valid when the function Q(x) has no zeroes in the interval of 
interest. Following Liouville and Green, we introduce a new independent variable 
t = f(x) implicitly determined as dt = \[Qdxjt. The choice is natural from the 
perturbation point of view, because even when du/d(x/e) is significant, dQ(x)/dx 
is not, so Q(x) can be regarded as a constant to order O(e ). Eq. (3.46) is thus 
converted to 

d 2 u „. ^du , A . 

w - u = 2eS ( X )-, (3.48) 

where S(x) = — (l/4:)Q~ 3 / 2 Q'(x) is assumed to be a slowly varying function on 



22 



the time scale t, of order unity, and S(x) 7^ for < x < 1. 

Naively expanding u as u(t) = uo(t)+eui(t)+- • •, we get the bare perturbation 
result 

t t 
u(t) =e*{ A + o4 y 5 - eA e- 2t ^ 5 (x(t')) e 2t ' dt'} 

to to 

t t 
+ e-^Bo + eB J S (x(t')) dt' - eB e 2t J S (x(t')) e~ 2 *' ' dt'} + 0(e 2 ), 

to to 

(3.49) 

where Aq, Bo are integration constants. The corresponding renormalized result 
is 



t 



t 



u{t) = e^Ai/j) + eAiji) J Sdt'} + e^js(^) + eB(fi) J Sdt'} + 0(e), (3.50) 



where O(e) refers to all regular terms of order e which remain finite even as 
t — to — > 00. The RG equation du/d/j = gives 



^ + e 1 -Q-^Q'(x^))C = 0(e% 



(3.51) 



where C = A or B. Again, equation (3.51) corresponds to the amplitude equation 
or slow motion equation. Setting \i = t and using dt = Q 1 / 2 dx/e, we get 



A(x) ~ Q- 1/4 (x), B(x) ~ OT 1 / 4 ^). 



(3.52) 



This is exactly the adiabatic invariant A(x)<5 1 / 4 (a:) = A(0)<5 1//4 (0) = constant. 
The physical-optics approximation for WKB equation (3.46) is recovered 



(x) -dQ-^ix) exp 



u 



as e — > 0. 



dx'^Qte 



+C 2 Q- lf Hx) exp 



dx'y/QW 

(3.53) 



The uniformly valid asymptotic result for the general linear boundary 
layer problem (3.44) is given by (3.45). For numerical evaluation of (3.33), we 
do not need any further expansion, because (3.45) is the uniformly valid result 
we want. To compare, however, with the conventional results due to asymptotic 
matching methods, let us make asymptotic expansions of Q(x). 
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As a simple check, let us assume that a(x), b(x) are some analytic functions, 
and a(x) > for < x < 1, and a(0) 7^ 0. Obviously, in the whole region 
< x < 1, as e — > 0, the term a 2 (x)/4e 2 is the dominant term, compared to 
a'(x)/2 and Simply Taylor expanding as 



, r— . a(x) e a'(x) b(x) 
2e 2 a(x) a(x) 



(3.54) 



and imposing boundary conditions y(0) = A, y(l) = £?, we obtain 



f'MSirff a(0) 
a(x) 



A -Be 



Jl a(O at - 



/•x 




a (€) 1 Mil 



(3.55) 



This expression can be simplified further, because the second term contributes 
appreciably only when x = 0(e 2 ) (e — > 0). Thus, 



y(x) ~ Sexp 



6(0 

«(0 



+ 



.4 -L'cxp / 



,-a(0) 2: /e 2 



(3.56) 



This is exactly the same as the uniformly valid leading boundary layer or WKB 
result. 

It is known that the case with a(0) = is subtle. For simplicity, we consider 
only the cases a(x) = x a , b(x) = 1, where a > — 1 so that there exists a boundary 
layer at x = 0. 

When a > 1, the thickness of boundary layer is of order 6 = O(e). When 
x ^> e, the term a 2 (x)/4e 2 dominates over other two terms, a'(x)/2 and b(x) 1 
in Q(x). However, when x ~ 0(e), we have to be careful with the asymptotic 
expansion of Q(x), since the dominant term now is b(x) = 1. Thus, as e — > 0, 
the leading term of >jQ(x) is 1. The final uniformly valid approximation is 



y(x) ~ B exp 



(l-x 1 -") 



a — 1 



+ A exp [— x/e] 



(3.57) 



When I a I < 1, it is straightforward to check that the boundary layer is of 
thickness of order 6 ~ e 2 ^ 1+a \ and that the first and second term in Q(x) are of 
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the same order, when x ~ 5(e). The uniformly valid expression turns out to be 



y(x) =Bexp 



a(0 



d£ 



+ 



A-Sexp / ^d£ 



X{ Q(0) ,-1/4 



exp 



X 

-/«{ 

. 



«(0 



a(0 , >/Q(0 



2e 2 



+ 



(3.5* 



Expanding the above leading uniformly valid result obtained with the aid of 
RG recovers the outer and inner solutions due to boundary layer theory and 
asymptotic matching. Note that the above results are obtained from the "inner 
expansion" alone without ever having to perform any asymptotic matching. This 
is practically important as we will see in the next section. However, the main 
message here is conceptually more important: conventional singular perturbation 
methods can be understood naturally as the standard renormalized perturbation 
procedure. 

3.7. WKB ANALYSIS II: TURNING POINTS 



In order to complete this section, we begin by presenting a general discussion 
of Schrodinger equations and one-turning-point WKB problem, and at the end 
of this section, we generalize the case to multiple-turning-point and multiple- 
boundary-layer problems. 

The Schrodinger equation which we will consider in this section is 



,d 2 u 
dx 2 



— Q(x)u(x), u(+oo) = 0. 



(3.59) 



When Q in (3.59) vanishes or changes its sign, the approach in the preceding 
subsection fails as can easily be seen from the presence of the factor Q' 1 / 4 . If Q 
has an isolated zero at x = of order a > 0, we can write locally Q(x) = x a ip(x) 
with a positive definite function ip without any loss of generality. A natural choice 
of the counterpart of the Liouville-Green transformation x — > t is to remove the 
zeros from dt/dx: we introduce a new independent variable t = f(x) implicitly 
determined as dt = ^/Q/t a dx/e, and integrating it gives 



2/(2+a) 



*(*) = 



2 + a 



dx'y/Q&j 



(3.60) 
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The original equation (3.46) is transformed into 



d u n „ , , s s du 



df , l"u ■ <S [!(,-))-. (3.61) 

where S = d[(t a /Q) 1/2 ]/ 'dx. Since t ~ x as x -> 0, 5 is a bounded function 
even near x = 0. Notice that in contrast to the conventional approaches due to 
Liouville and Green or Langer, we do not introduce the transformation for the 
dependent variable, which will be produced by the RG procedure. Here we work 
out the simplest case a = 1. 

Expanding naively u in powers of e as u = uo + eu\ + e 2 U2 + • • •, we obtain 
the bare perturbation result to order e, 

t t 
u = C Ai(t) - eC 7r{Ai(t) J dt'S{t')Ai'(t')Bi(t') - Bi(£) J dt' S{t')M{t')A]!{t')}, 

to to 

(3.62) 

where Ai, Bi are two linearly independent Airy functions, and the Bi(t) function 
is already discarded in the zeroth order solution, since it does not satisfy the 
physical condition u(+oo) = 0. In the limit t — to — > +oo, the second term of 
the first order perturbation Bi(t) dt' S(t')Ai(t')Ai' \t') remains finite. However, 

the term J to dt'S(t')Ai(t')Bi(t') diverges and must be renormalized, giving the 
renormalized perturbation series 



u 



= Ai(£) | C(/jl) - eC(ii)iv J dt'S(t')Ai'(t')Bi(t') J + 0(e), (3.63) 



where C((i) is the counterpart of the bare amplitude Co(to) 5 an d 0(e) refers to 
all finite regular terms of order e even in the limit t — to —> oo. The RG equation 
du/dfi = gives 

^j^- + eC(/j)nS(ii)Ai'(ii)Mv) = 0(e 2 ). (3.64) 
Integrating (3.64) and setting \i = t, we get 

C(t) = O(0) exp | - tt j dt'AS'{t')Bi{t')^{\og[{t'/Q) 1 l 2 ] }|, (3.65) 

where Oq is a constant of integration to be determined by boundary condition at 
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t — 0. Thus we have arrived at the adiabatic invariant 



C(t) exp /tt J dt'ki!{t')Bi{t')^[\og [it'/Q) 1 ' 2 ] } j, (3.66) 



which differs from that usually obtained, leading to the the final uniformly valid 
solution 



t 



u 



= C(0)exp{-7r J dt>Ai\t>)Bi(t>)-^(\og[(t'/Q) 1/2 ])}Ai(t), (3.67) 



o 

where t(x) = /* dx'y/Q(x'f) . 

The RG result (3.65) differs from the standard Langer formula, since (3.65) 
involves Airy functions Ai and Bi. Note that the new variable t given in (3.60) 
is a function of e, and that as e — > for fixed x, and t — > oo. In this limit, we 
can resort to the asymptotic properties of the Airy functions Ai(i) and Bi(t) for 
t — > oo, and find that Ai'(£)Bi(f) ~ — l/2n, as t — > oo. Thus, (3.65) recovers the 
standard result 

C(t(x)) = C(0)(t/Q) 1 /\ (3.68) 

However, the RG equation (3.64) is valid not only for relatively large fi, but 
also for small \i. For this reason, we expect that (3.67) is a better uniformly 
valid approximant than the standard Langer formula, for small and intermediate 
values of t, or for relatively large (or not small) e cases. This is verified and can 
be clearly seen in Fig. 1, where we compare the RG result (3.67), the standard 
Langer formula, and the exact numerical solution of equation (3.46) for several 
values of e. Thus, the RG results (3.67) without asymptotic matching improve 
upon those obtained by the standard analysis. 

To conclude this section, we briefly outline the recipe to generalize the meth- 
ods for multiple-turning-point and linear multiple-boundary-layer problems. (For 
linear cases, with the help of the transformation (3.45) both problems can be 
transformed into the canonical form and can be treated in a unified way.) We 
need only consider the case in which Q(x) in (3.46) has multiple turning points. 
Without loss of generality, we assume Q has the form: Q(x) = f(x)^(x), where 
f(x) = (x — x\){x — X2) • • • (x — x n ), n > 1 is a polynomial of x with n zeros 
x\ < X2 < ■ ■ ■ < x m and ip(x) > has no zeros. The general strategy is first to 
introduce a new independent variable t defined implicitly as dtj y/Q/ f(t)dx/e, 
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where / is chosen to cancel all the zeros of Q. Then we develop the straightfor- 
ward perturbation series for the resultant equation, and renormalize the integra- 
tion constant to absorb the secular divergence. This procedure avoids performing 
multiple connection formulae matching and leads to a uniformly valid approxi- 
mation. For higher order WKB problems or linear boundary layer problems, the 
generalization of the methods given here is straightforward. 



4. Switchback Problems 



In previous sections, we have already seen that the RG approach not only 
has conceptual, but also technical advantages compared with various conventional 
methods. In this section, we will demonstrate this further, by studying, with the 
aid of RG more complicated problems which involve the so-called 'switchback'. 
In switchback problems, as conventionally treated, only through subtle analysis 
in the course of actually solving the problem is it possible to realize the need for, 
e.g., unexpected order terms to make asymptotic matching consistent. 



4.1. Example 1: Stokes-Oseen Caricature 



A model example is a caricature of the Stokes-Oseen singular boundary layer 
problem, which describes the low Reynolds number viscous flow past a sphere of 
unit radius. The main result of this problem has been presented in ref. 16, and 
in the following, we will only briefly summarize the final results and make some 
additional comments. 

5 

The equation is 

d 2 u 2 du du . . . , . A . 

*T + 7* +e "*-°- «(l)=0,«(oo) = l, (4.D 

where e, the Reynolds number, is a small non-negative constant. This is a very 
delicate singular boundary layer problem, with complicated asymptotic expan- 
sions and matching, involving unexpected orders such as elog(l/e). 

Since there exists a boundary layer of thickness 5 = 0(e) near r = oo, setting 
x = er transforms (4.1) into the following 'inner' equation: 

d 2 u 2 du du , . n . . . , _ 

— T H — + u— = 0, u[x = e) = 0,u(x = oo) = 1. (4.2) 

dx z x dx dx 

Using RG theory, the final uniformly valid result is found to be, to order 
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Ai = l/e2(e), 

u(r- e) = 1 - e 2 (er)/e 2 (e) + o{ [l/e 2 (e)] 2 }, (4.3) 

where the exponential integral e 2 (t) = J t dpp~ 2 e~ p , whose asymptotic expansion 
as t — > is given by e 2 (t) ~ l/£+log£+ (7-I) -t/2 + 0(t 2 ) with Euler's constant 
7-0.577---. 

5 

The result from asymptotic matching is given by the following expression. 
For r fixed, we have, to O (e 2 log 2 (l/e)) , 

u(r) ~ (1 - 1) + elog(l/e)(l - i) + e[ - logr + (1 - 7 )(1 - 1 - ^)], (4.4) 

while for p = er fixed, to O (e 2 log(l/e)), 

u(p) ~ 1 -ee 2 (p). (4.5) 

Accordingly, examining the asymptotic result of (4.3) in the limit e — > 0, by 
expanding both e 2 (er) and e 2 (e) for r fixed, and e 2 (e) only for p = er fixed, re- 
spectively, it is found that the resulting asymptotic solution using RG is correct to 
0[e log(l/e)] and agrees with that obtained by asymptotic matching. Note that 
in our method, the eloge term appears naturally from the asymptotic expansion 
of e 2 (e), whereas some artistry is required to obtain this term conventionally. To 
recover the 0[e] term with (log r)/r, we have to extend the RG calculation to 
order 0[(l/e 2 (e)) 2 ]. Thus, the result to 0[e] given by asymptotic matching 5 is 
obtained from the renormalized perturbation expansion to 0[(l/e 2 (e)) 2 ]. 

This fact may suggest that our RG result is inferior to the conventional 
one. It is important to notice, however, that neither the asymptotic expansion 
(4.4) augmented with the (log r)/r term of order e nor (4.5) is uniformly valid 
in its variable r or p, respectively. In contrast, it seems that our full result 
1 — e 2 (er)/e 2 (e) to order Ai = l/e 2 (e) is uniformly valid as is clearly seen in Fig. 
2. 

As discussed in the preceding paragraph, (4.3) is not an asymptotic series 
in powers of e; thus, one might conclude that our result is not even an asymp- 
totic series in any sense. Recall, however, that the asymptotic expansion of a 
function is unique only when an asymptotic sequence of functions is fixed. The 
choice of the sequence is a question of vital importance, if one wishes to have 
a useful asymptotic series. In the conventional singular perturbation methods, 
an asymptotic sequence is selected by the matching conditions. However, there 
is no compelling reason to believe that the selected sequence is practically the 
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best asymptotic sequence (of course, it should be the most convenient one for 
the matching procedure). As we have seen, the RG approach also produces an 
asymptotic sequence {Aj(e)} from the requirement to satisfy the boundary con- 
dition order by order. Therefore, we propose the point of view that a consistent 
and presumably better asymptotic expansion (starting with Ai = l/e2(e) in the 
present problem) may be obtained by RG. The standard e expansion may well be 
an inferior asymptotic expansion to our expansion. In addition, the superiority 
to the RG approach can also be seen from the fact that a closed expression uni- 
formly valid for the whole (infinite) interval has been obtained for the problem, 
which is not the case for the standard asymptotic matching method. 

4.2. Example 2: Difficulty with Asymptotic Matching 

To illustrate that the RG method is generally simpler to use, and yields prac- 
tically better approximants than other methods, let us next consider a 'terrible' 
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problem whose model equation can be written as 



where e is a small non- negative constant, and a = or 1. For a = 1, the 
asymptotic matching is notoriously difficult, because an infinite number of terms 
must be calculated before even the leading order can be matched successfully. 
We will see how the RG avoids such difficulties in obtaining the leading order 
result uniformly valid for the entire interval 1 < r < oo. 

Since there exists a boundary layer of thickness 5 = 0(e) near r = oo, setting 
x = er transforms (4.6) into the following 'inner' equation: 



As in other boundary layer problems, let us first look for the general form of 
the solution, and then impose the required boundary conditions to determine the 
constants of integration left in the solution. To do so, we solve (4.2) as an initial- 
value problem, given an initial condition u(xo) — Aq at some arbitrary point x = 
xq, where Aq is a finite constant. Assuming a naive expansion u(x; e) = uq(x) + 
Xi(e)ui(x) + A2(e)«2(x) + ■ • • with initial conditions uo(xo) = Aq, Ui(xo) = 0, i = 
1, 2, • • •, where the asymptotic sequence Aj(e), i = 1, 2, ■ • • are to be determined 





(4.6) 




u(x = e) = 0, u(x = oo) = 1, (4.7) 
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later, we obtain 

d 2 uo 1 duo ( duo \ ^ duo 



dx 2 ^ x dx ^ \ dx J ^ dx ^ ^ 

The finite uniform solution can be guessed as uq(x) = Aq, because the uniform 
field should not be affected appreciably by the distant disturbance source. Thus, 
the goal is to find out the small perturbation effect on this uniform field in the 
presence of a distant disturbance. 

The equation for u\ is 

d 2 ui ( 1 . \ du\ n . . _ 

+ " + ^0 -rr = 0. (4.9) 



dx 2 \x ) dx 

We easily see that the equation satisfied by U2 which is significantly different 
from (4.9) (i.e., with a forcing term) appears only if \ 2 /\2 = 0(1). We will show 
that indeed the choice A2 = X 2 works. The nontrivial equation at order A2 = X 2 
can be written as 

d 2 U2 fl . \ du2 ,dui*o du\ /A in . 

z + [- + A — ± = -a(— -f-ux— -. (4.10) 



dx 2 \x J dx dx dx 

The perturbation result is given by 



u(x) =A + X(e)A w [e 1 (Aox ) - ei (A x)] + X 2 (e)^A 2 [ei(A x ) - e^Aox)] 

- a-Al [ei(A xo) - ei(A x)] 2 - A\qA q 1 [eo(A x)ei(Aox) 

- 2ei(2A x) - ei(A xo)eo(Aox) + eo(Aoxo)ei(A x) 

- e (Aoxo)e 1 (A xo) + 2ei(2A x )] } + 0[A 3 (e)], 

(4-11) 

where the exponential integral e\(t) = f t p 1 e p , Ai(e) is already replaced by 
A(e), and A\o, A2 are constants of integration. When xo is very small and x — xq 
is large, the divergence arises from those terms containing e\(Aox) or e\(AoXo), 
but not eo(Aox) or eo(AoXQ). To remove the divergence from these cross terms of 
eo and e\, presumably both Aq and A\ must be renormalized. The renormalized 
perturbation result reads 

u(x) =A(u) + X(e)A 1 (u)[e 1 (Au) - ei(Ax)] + X 2 (e)^A 2 [e 1 (Au) - ei(Ax)] 

- a-A\[ei(An) - ei(Ax)f - A\A~ x \eo(Ax)(e\(Ax) - ei(Au)) 

- 2e x (2Ax) + 2 ei (2Au)] } + 0[A 3 (e)], 

(4.12) 
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where A(n), Ai(/j,) are finite counterparts of Aq(xo), Aiq(xq), and \x is some arbi- 
trary length scale. The RG equation du/dfi — gives 

^ = -\(e)aA\n- l e- A » + 0[A 2 (e)], (4.13) 

^ = -X(e)A 1 ^- 1 e- A ^ - \ 2 {e)A\A- 1 ii- l e- 2A ^ + \ 2 {e)A 2 ^ 1 e^ + 0[A 3 (e)]. 

(4.14) 

Now we discuss the a = and a = 1 cases separately. For a = (4.13) 
suggests that A\ can be treated as a constant and there is no need to renormalize 
it. Solving (4.14) to order A(e) and setting \x = x and x = er in (4.12) , we obtain 

u(r) = 1 - X(e)A 1 ei(er) + A 2 (e), (4.15) 

where use is already made of the boundary condition u(r = oo) = 1. Impos- 
ing u(r = 1) = determines X(e)Ai = l/ei(e) from which A(e) can be cho- 
sen as A(e) = l/ei(e), whose asymptotic expansion in the limit e — > 0+, is 

A(e) ~ l/log(l/e) + 7/log 2 (l/e) H , giving all necessary orders required in the 

asymptotic matching. Accordingly, A\ = 1. Thus, the uniformly valid asymp- 
totic result can be written in a single expression as 

u(r) ~ 1 - ei(er)/ei(e) + 0[(l/ ei (e)) 2 ]. (4.16) 
For a = 1 solving (4.13) and (4.14) to order A(e), we get 

AM = : xm4 ( ?\ , S + 0[>?(e)}, ( 4 - 17 ) 
1 - A(e)Ai(oo)ei(//) 

A(^) = log {l - A(e)A 1 (oo)ei(^)} + A(oo) + 0[A 2 (e)], (4.18) 

where Ai(oo), A(oo) are constants of integration to be determined by the bound- 
ary conditions. Setting [i = x and x = er in (4.12) we have 

u{r) = log {l - A(e)Ai(oo)ei(er)} + A{oo) + 0[A 2 (e)]. (4.19) 

Using boundary conditions u(r = oo) = 1 and u(r = 1) = produces A(oo) = 1 
and A(e)Ai(oo) = (1 — l/e)/ei(e). Again we may choose A(e) = l/ei(e), and then 
Ai(oo) = 1 — 1/e. Finally the uniformly valid asymptotic result is given by 

u(r) ~ log [l + (e - l)ei(er)/ei(e)j + 0[(l/ei(e)) 2 ]. (4.20) 

Comparing the RG results (4.16) and (4.20) and the corresponding asymp- 
totic matching results, again we find the RG results are more accurate. 
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5. Reductive Perturbation Theory and RG 



In previous examples, we have already mentioned the idea that amplitude 
or phase equations are RG equations. We will demonstrate that the RG theory 
is a general and systematic method to derive slow motion equations, even for 
those complicated problems for which no explicit analytic zeroth-order solutions 
are known. In previous reports we already discussed the one-dimensional Swift- 
Hohenberg equation 16 and the Burgers equation 11 as renormalization group equa- 
tions. Center manifold theory can be considered from the reductive perturbation 
point of view, because it also extracts slow motion equations on the manifold. 
Thus, we may expect that the center manifold theory can also be interpreted as 
an application of the renormalization approach as well. 



5.1. Newell- Whitehead equation 



The example we consider here is the two-dimensional Swift-Hohenberg equa- 
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tion widely used as a simple model of the Rayleigh-Benard convection, 



2 \ 2 

«*"« 3 - (^2+1^2+^ «, (5-1) 



du 3 / d 2 d 

dt \dx 2 dy 



where e is a control parameter or a reduced Rayleigh number, a measure of the 
degree of convective instability of the stationary state u = 0. For small positive 
e, the system exhibits a supercritical bifurcation. Since we wish to treat eu — n 3 
as a perturbative term, to be consistent eu and it 3 must be of the same order. We 
scale u as -Jeu, and denote the new u with the same symbol. Then, the original 
equation reads 

9« ■ .3, ( a 2 . a 2 , ,S 1 



W = £(M - u) -U? + 5? + *J (5 ' 2) 

We consider this in the whole plane for all positive t. As a zeroth order solution, 
we choose the roll solution along the y-axis: Ae tkx + complex conjugate, where 
A is a complex numerical constant. We expand u around this solution as u = 
Ae %kx + eu\ + h complex conjugate. The first order correction obeys 

^ + g + 5 + ^) 2 « 1 = ( 1 -3 W Ve- (5.3) 
Here, to study only the singular behavior of u\, e 3tkx and similar non- resonant 
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terms are ignored. We rewrite this equation as 

[Li + L 2 + L 3 + L 4 ]ui = (1 - 3|A| 2 )Ae fc , (5.4) 
where the operators are given defined as 

£ ' s «- i ^(a? + fc J • ^H^W £ ^v' (5 ' 6) 

We must look for space-time secular terms in the solution. Secular terms appear 
only in the special solution of the equation consistent with the inhomogeneous 
term. In order to find (space-time secular) special solutions of (5.4) we have only 
to solve LiUsi = (1 — 3|A| 2 )Ae* fc:r separately, and to make the linear combination 
of their solutions as ^^usi with = 1. This is because all four operators 
Li commute, and Lje lkx = 0, so that LiLjusi = 0. A trivial calculation gives 

usi = tA(l-3\A\ 2 )e ikx . (5.6) 

us2 is governed by 

<}1 + k 2 ) 2 u S 2 =0- + ik] 2 0- - ik] 2 u S2 = (1 - 3\A\ 2 )Ae^. (5.7) 



dx 2 J \dx J \dx 

That is, 



{r x -' k ) »«=->-viv b (5.8) 

Here we do not pay attention to inhomogeneous terms nonresonant with the 
operator. Hence, the most singular part is 



u S 2 = -^A(l-3\A\ 2 )e^. (5.9) 



Similarly, we get 



and 



u S3 = ^-A(l-3\A\ 2 )e ikx : (5.10) 

u SA = - 3\A\ 2 )e ikx . (5.11) 

In this way we get the following perturbation result, 

u = Ae ikx +e - + + pJL^j A(l-3\A\ 2 y kx + c.c+- ■ ■ . (5.12) 

Here all the less singular terms (e tkx times 1, x, y, xy, y 2 , y 3 ), higher order terms 
and nonsecular terms (those terms which do not grow indefinitely far away or in 
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the long future) are omitted. These terms will not contribute to the final result, 
as shown in the argument below. Now, the secular terms are absorbed into the 
redefinition of the amplitude A as follows. We introduce regularization points 
X, Y and T and split, for example, x a as x a — X a + X a (for some exponent a), 
and absorb X a into A. Thus we get, 

u =A(X, Y, TY k * + e L(t - T) - ^ ~ k P + ^ ~ z f ^ 

x (5-13) 
+ /"4^^ A(l - 3| A\ 2 )e ikx + ■■■. 

Since u should not depend on X, Y or T, the renormalization group equation, 
to 0(e), reads d a+l3+1 u/ 8T a dXl 3 dY 1 = for any positive integers a, f3, 7 with 
aj3^f 7^ 0, where values of a, ft, 7 are chosen in such a way that the universal slow 
motion equation we are seeking is independent of any system details. Thus, we 
have 

OA 

__ emA (l-3|A| 2 ) = 0, 
f) 2 A 1 

" +^ 2 _4(1-3|A| 2 ) = 0, 



dX 2 ' ^4k 2 ' 
d 3 A 1 . ... 2 



(5.14) 



oxdY 2 ™uk 
d 4 A 



e ^Hk A{1 ~ m ) = ' 

- e/i 4 A(l - S\A\ 2 ) = 0. 



dY 4 

Obviously, m are still almost arbitrary and must be fixed by the auxiliary con- 
ditions. Therefore, to get an auxiliary condition free equation of motion, we use 
^2 H-i = 1 to arrive at the following RG equation after equating X, Y, T and x, y, t, 
respectively 

6 4 + {~^&2 + «^ + * = ^(1 - 3|^| 2 ). (5.15) 

Thus, we have arrived at the Newell- Whitehead equation. 

Let us compare this derivation with the conventional method, for which a 
summary may be found in the Appendix to the review article by Cross and 
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Hohenberg. Perhaps the most notable point is that no scaling of spatial vari- 
ables like x — > e 1//2 x, y — > e 1//4 y is needed. Furthermore, the expansion is a 
straightforward one in terms of e instead of e 1 / 2 . That is, the result is almost au- 
tomatically obtained from the global well-definedness of the perturbation result. 
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If there are no spatial degrees of freedom, each step of the standard reduc- 

7 

tive perturbation using the solvability condition and that in the RG derivation 
above are in one-to-one correspondence. However, if there are spatial degrees of 
freedom, the standard reductive perturbation regards the spatial derivatives as a 
perturbation if the zeroth order solution is space-independent, or uses the multi- 
ple scale analysis if the zeroth order solution is spatially varying. In contrast, in 
our RG approach, spatial and time coordinates are treated on an equal footing, 
and the correct scalings of variables are given automatically. 

As the reader may have realized, kinetic equations are expected to be deriv- 
able as slow motion equations from the BBGKY hierarchy. For example, the 
Boltzmann equation can be derived by an RG method. Thus we suggest that it 
is a rule that slow motion equations are RG equations. 



5.2. Center Manifold and RG 

In this section, we discuss briefly the general relationship between RG theory 
and center manifold theory. In the general theory of reduction, we wish to know 
the slow manifold (e.g., inertial manifold, center manifold) which attracts all the 
long-time asymptotic solutions, and the equation of motion on the manifold. It 
is well known that the center manifold reduction and normal form theory 6 have 
played a significant role in studying instabilities and bifurcations encountered in 
dynamical systems and fluid dynamics. In many circumstances, this approach 
provides a greatly simplified picture of complicated dynamics by reducing the 
dimension of the system without losing essential information concerning the in- 
stability and bifurcation. In addition, the local dynamics on the center manifold 
constructed in this way is invariant or universal, in the sense that the structure 
of the reduced system is independent of specific physical models under considera- 
tion. Thus a variety of different phenomena can have the same type of bifurcation, 
belonging to the same universality class in the parlance of RG. Although the cen- 
ter manifold fits in the RG picture clearly, the general correspondence between 
them has not yet been established. In certain cases such as the weakly nonlinear 
stability of fluid motion, the equivalence of the method of center manifold, the 
method of multiple scales, and the method of amplitude expansion has been es- 
tablished explicitly by applying these methods to the derivation of the Landau 
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equation from the Navier-Stokes equation to the seventh order. 
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To illustrate the relevance of RG, let us consider the following set of equations: 



dx , . 

37 =f(x,y), 

f (5.16) 

^ = -y + </(*,y), 

where / and g are higher order in the sense that f(Xx, Xy) or g(Xx, Xy) is 0(X 2 ) 
for small A. Thus the variable y decays quickly but x does not. Hence, the long 
time behavior of the system is expected to be confined close to a local 1-manifold 
near the origin. This local manifold is the center manifold (not unique) , and the 
long time behavior of the system is governed by the equation of motion defined 
on this manifold. Thus, as discussed at the beginning of this section, the problem 
of finding a center manifold and the equation on it is a problem of extracting 
slow motion behavior of the system. In this sense, this problem and the general 
reductive perturbation can be treated in a unified fashion. Since we are interested 
in the local center manifold, we may rescale the variables as x — > Xx and y — > Ay, 
and may assume that A is small. Therefore, instead of the original system (5.16), 
we study 

f =Xf(x,y), 

dy , . 

- = - y + Xg(x,y). 

We assume the following formal expansions: 



(5.17) 



f{x, y) =hox 2 + fuxy + f 02 y 2 + X(fs x 3 + f 2 \x 2 y + fi 2 xy 2 + fay 3 ) H 

g(x, y) =92ox 2 + gnxy + gmy 2 + A(#30^ 3 + 92ix 2 y + guxy 2 + gosy 3 ) H • 

(5.18) 

The standard approach goes as follows: Let y = h(x) be the formula for a center 
manifold. Then we get the following differential equation for h: 

-h(x) + Xg(x, h(x)) = Xh'(x)f(x, h{x)). (5.19) 

This equation is usually solved by perturbation: h{x) = Xh 2 x 2 + X 2 h^x 3 + • • •. 
The result is 

V = ><92ox 2 + X 2 [g 20 (g n - 2/ 20 ) + Sso]^ 3 + 0[X 3 ]. (5.20) 

The equation of motion on the center manifold is obtained by substituting y with 
h(x) in the equation for dx/dt. 
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Our RG program starts with the construction of a power series expansion 
of the solution for (5.17) in terms of A as x = xq + Xx\ + X 2 X2 + • • •, and y = 
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Vo + Xyi + X y 2 + ■ ■ ■■ Paquette has also pursued the same line independently. 
A lengthy but straightforward calculation gives 

x =A + Xf 20 A 2 t + X 2 (f 2 A 3 t 2 + fn92oA 3 t + f 30 A 3 t) 
+ X 3 {t 3 f 3 A 4 + !(/ii/ 20 <720 + hof2o)A 4 t 2 

+ [-2^20/20/11 + hmogn + hmo + /2o^lo + /2i^20 + Uo]A 4 t ( 5 - 21 ) 
+ CT}+0(A 4 ) 
y =Xg 2Q A 2 + X 2 [2g 20 f 20 A 3 (t - 1) + ^n^ 3 + 9wA% 

where CT denotes the constant terms and A is the initial condition for x. Here 
we have discarded all the exponentially decaying terms. For example, to the first 
order the full solution reads 

xi = f 20 A 2 t - fuAoBoe-t - \h 2 B 2 e~ 2t + A u 

yi = 92oA 2 + gnAoBote- 1 - go 2 B 2 e~ 2t + B x e~\ 

where Aq, Bq, B\ are numerical constants dependent on the initial data. The 
exponentially decaying terms do not contribute to the secular behavior of per- 
turbation series. We absorb the secular terms proportional to the powers of T 
into the redefined A by splitting t as t a — T a + T a , where a is an appropriate 
integer. That this can be achieved consistently must be checked order by order. 
The simplest way may be to introduce the renormalized counterpart Ar of A as 
A = Ar(1 + Xui + X 2 uj 2 + X 3 ujs + - • •), where u>i are determined to remove the pow- 
ers of T from the perturbation result for x after splitting t. The renormalization 
condition can be written as 



(5.23) 



A R (1 + XuJi + X 2 uj 2 + X 3 u 3 + •••)+ A/ 20 (l + Xui + X 2 u 2 ) 2 A 2 t 
+A 2 (/ 2 V 2 + hmot + hot)A 3 (l + XiV!) 3 + A 3 [- • •] = A R . 
From this, order by order in powers of A, we can fix Ui as 
u>i = - f 20 At, 

u 2 =fl A 2 t 2 - f n g 20 A 2 t - f 30 A 2 t, 

us =A 3 | - /| t 3 + ^(/n/20020 + /3o/ 2 o)t 2 (5.24) 

+ [2y20/20/ll< - /llfi r 20fi r ll - /ll#30 - /20^20 ~ ^lfl^O - f40}t 
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The renormalization group equation reads 



^ = ^{Ar(T)(1 + Xcv, + X 2 u> 2 + A 3 W3 + •••)} = 0. (5.25) 

Introducing the explicit forms of uii into this equation, we experience almost 
miraculous cancellations of all the terms containing powers of t explicitly to have 

^ =XA 2 R f 20 + X 2 A%(f imo + ho) (526) 
+X s A 4 R (f 11 gso + f 2 ogi + hig-xs + ho - ho f 11920) + C[A 4 ], 

where t is identified with T. This agrees with the conventional result. For y, 
after renormalization, all the explicitly t dependent terms disappear to order A 2 , 
and 

y = Xg 20 A 2 R + X 2 (g 20 gn + g 3 o - 2f 20 g 20 )A R + 0[A 3 ]. (5.27) 

This also agrees with the result given above. 

The formal solution (5.21) is order by order in A obtained from the true so- 
lution by discarding the transcendentally small terms in the large t limit. Notice 
that in x n the highest power of t is n (for y n it is less), so that up to a given 
order n, by choosing A such that At = 1, we can make the contribution of the sum 
of the transcendental terms (such as e _1//A ) less than any small positive number 
for sufficiently large t. In this way, locally up to any finite order in A, the series 
obtained as the singular (or non-decaying) terms describes the asymptotic be- 
havior of the system. Therefore, if the system has a unique solution to the initial 
value problem (near the origin), then we can uniquely determine these series, 
and they give a parametric representation of an approximate center manifold. In 
the present context, renormalizability means that the motion on the approximate 
center manifold is autonomous. The renormalization reorganizes the expansion 
so that dx/dt is not explicitly time-dependent. 

The RG procedure given above is actually much more tedious than the con- 
ventional approach. However, the obtained center manifold by RG need not be 
expandable in terms of x. Thus, the RG method works in some cases even when 
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the conventional approach is not applicable. 



39 



6. Summary 



In this paper, we have demonstrated that various singular perturbation meth- 
ods and reductive perturbation methods may be understood in a unified fashion 
from the renormalization group point of view. Amplitude equations and phase 
equations describing the slow motion dynamics in nonequilibrium phenomena 
are RG equations. The RG method seems to be more efficient and simpler to 
use than standard methods. This is practically meaningful, because RG could 
yield superior approximations without using often tedious asymptotic matching 
techniques. 

Probably the most outstanding question is to justify mathematically the gen- 
eral renormalized perturbation approach developed in this paper. The rigor- 
ous and constructive renormalization group approaches of Bricmont and Kupi- 
ainen and our formal perturbative approaches have almost no common technical 
ground, although their philosophy is identical. Consequently, we do not have 
even a hint as to how to rigorize, or estimate the errors of our approach. 

The Wilson-style RG 44 and Bricmont and Kupiainen's related constructive 

renormalization group approaches 13 can be implemented numerically, 45 especially 
with the interpolation-resampling scheme which produces a 'virtual continuum' 

46 

to allow smooth scaling of any function on discrete grids. 
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FIGURE CAPTIONS 

Fig. 1. Comparison of the RG result (3.67), the standard Langer formula, and the 
numerical solution of equation (3.46) for e = 0.5 and e = 1.0. 

Fig. 2. Comparison between the numerical solution of eq. (4.2) for several values 
of e, the first order RG result 1 — e2(er)/e2(e), and two matched asymptotic 
expansions (one at fixed r, the other at fixed p = re), as derived in ref. 5. 
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